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Abstract 

We analyze the reheating in the modification of the vMSM (Standard Model with three right handed neutrinos with masses 
below the electroweak scale) where one of the sterile neutrinos, which provides the Dark Matter, is generated in decays of the 
additional inflaton field. We deduce that due to rather inefficient transfer of energy from the inflaton to the Standard Model sector 
reheating tends to occur at very low temperature, thus providing strict bounds on the coupling between the inflaton and the Higgs 
particles. This in turn translates to the bound on the inflaton mass, which appears to be very light 0.1 GeV < mi < 10 GeV, 
or slightly heavier then two Higgs masses 300 GeV < mi < 1000 GeV. 
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1. Introduction 

In |H Q] it was shown that within the Standard Model 
(SM) complimented with three right-handed neutrinos 
Ni with the masses smaller than the electroweak scale 
one can simultaneously explain both the dark matter 
and the baryon asymmetry of the universe 
SBSiSElHEi-This model dubbed as ^MSM 
represents a particular realization of the seesaw exten- 
sion of the SM and is fully consistent with the current 
experimental data from the light neutrino sector. How- 
ever, generation of the proper Dark Matter abundance 
of the sterile neutrino is not simple during the thermal 
evolution of the Universe, and requires some amount of 
fine-tuning lfl3l [fl . Being very weakly coupled, sterile 
neutrinos do not reach thermal equilibrium, so an inter- 
esting possibility is to generate them before the begin- 
ning of the standard thermal history. In |@] such mech- 
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anism was proposed, where the i^MSM model was ex- 
tended by adding the inflaton field, which generates all 
the masses in the model and decays into the SM parti- 
cles and sterile neutrinos after inflation, 



fi 

£„MSM (£*/MSM[Mi->-0] - -TtNjNtX + h.C.) + 



:(^x) 2 + l^$| 2 -y(<f>,x) , 



(1) 



where $ and \ are the Higgs and the inflaton fields 
correspondingly and 



CvMSM[Mi^o] = £msm + Niidf/y^Ni 

- F aI L a Nr<S> + h.c. 



(2) 



is the ^MSM Lagrangian with all the dimensional pa- 
rameters being put to zero. The potential V(Q, x) i a 1 1 



1 In order to avoid the domain wall problem a cubic term p,\ 3 can 
be introduced. It will be further assumed that p < ^/a 3 /A i'ew ■ In 
that case such term has no influence on the dynamics of the model 
during the reheating stage, and the relation j4) for the values of the 
parameters considered in the Letter is not altered significantly either. 
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i TO 2 x 2 + y , (3) 



where Vq = -jS- was introduced in order to cancel the 
vacuum energy. Expanding (0) around its vacuum ex- 
pectation value one has the relation between the infla- 
torOJmass mi and the Higgs mass mn- 



mi 



(4) 



If a > (3/2 the inflaton mass is smaller then the Higgs 
mass and, therefore, the decay of the inflaton into the 
Higgs can only occur in a thermal bath. In what follows 
we will first concentrate on this particular case. Param- 
eter (3 is fixed by the COBE normalization of the ampli- 
tude of scalar perturbations IU5I1 . (3 ~ 1.3 x 10 -13 . Pure 
quartic potential inflation is currently disfavored by the 
WMAP5 data 0] because of the too large predicted 
value of the tensor to scalar amplitudes ratio. However, 
if one allows non-minimal coupling of the inflaton to 
gravity Il7tl one can bring this potential in agreement 
with the data. This, in turn, will influence the bounds 
on the inflaton mass. We will discuss this in the end of 
the Letter. 

The upper constraint on the value of a comes from 
the requirement that radiative corrections do not spoil 
the flatness of the inflaton potential and is given by a < 
3 x 10~ 7 . This corresponds to the lower bound on the 
inflaton mass 



mi > 0.07 



m H 



150 GcV 



1.3 x 10- 13 



GeV . (5) 



One should note that larger values of a (leading to 
smaller inflaton masses) may also be possible, but then 
the analysis of the loop corrections to the effective po- 
tential of the inflaton becomes important. 

The lower bound on a comes from the requirement to 
have successful baryogenesis in ^MSM [2|. To allow for 
efficient sphaleron conversion of the lepton asymmetry 
to baryon asymmetry requires the reheating temperature 
to be larger then roughly 150 GeV JH]. In |Q| it was 
advocated that the resulting lower bound is a > (3 ~ 
10~ 13 . Below we will argue that the lower bound is 
quite a bit stronger which leads to a narrow window for 
the inflaton mass. 



2 Notations / and H will be used throughout the Letter to represent 
the diagonalized excitations above the vacuum expectation value for 
(3)- / is the one mostly mixed with inflaton \, and H mostly mixed 
with the SM Higgs <£. 



2. Reheating bounds 

Reheating after inflation proceeds through a regime 
of the parametric resonance. The dynamics of the mod- 
els with potentials similar to (O in the parametric reso- 
nance regime was studied via analytic methods in, e.g. 



H19L 12011 . The analysis of the late stages of preheating 
was made possible with the lattice simulations package 
LatticeEasy lr2ll I22I I23I l24ll . In particular, the preheat- 
ing in the model with the potential which contains only 
first two terms in (01 have been studied in 12211 . 

At large values of the inflaton field \ the behavior is 
that of the pure quartic inflation. The expectation value 
of the Higgs field $ is set along the flat direction: | $ | 2 — 
jX 2 . After the end of inflationary slow roll regime the 
inflaton field starts to oscillate. In the very beginning all 
the energy is stored in the zero (or homogeneous) mode 
of the inflaton \o, an d a H other modes are absent. The 
oscillations of \o initially excites the nonzero modes 
of both the Higgs and the inflaton. One can compare 
the contribution of the zero mode of the inflaton to the 
effective masses of the Higgs and the inflaton: 



off,* 



l eS, X 



(6) 



If a > (3 the corresponding contribution to the effective 
mass of the Higgs is larger. Therefore at early stages 
of the evolution the energy transfer into the Higgs par- 
ticles is the dominating process. This is in accord with 



II19LI20I1 . and can be inferred from the early time behav- 
ior of the number densities shown in Fig.Q] One could 
then expect that the whole energy of the inflaton field 
will be transferred exponentially fast to the Higgs parti- 
cles[f] Since the Higgs decay to the SM fields and their 
consequent thermalization are fairly fast compared to 
the Hubble rate one could then estimate the resulting 
reheating temperature as in |0] 



Tr ~ mpi 



a 

_g*A 



(7) 



which for A ~ 0.1, the number of the SM d.o.f. 3* ~ 
10 2 and a > (3 leads to the values of Tr which greatly 
exceed the freeze-out temperature of the sphaleron pro- 



3 One can easily verify that only a small fraction of the energy 
of the inflaton is drawn into sterile neutrinos because of the small- 
ness of the Yukawa couplings fj (at most ~ f0~ 7 for the heaviest 
sterile neutrinos lighter, then the inflaton). In particular, the rate 
r(J — > NN) typically equilibrates at the temperatures below the 
temperature of the electroweak phase transition. The process involv- 
ing the SM Yukawa couplings F a j proceeds via the Higgs particle, 
and is even more suppressed. 
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Fig. 1. Number densities of the Higgs boson and the inflaton are 
shown for different values inflaton-Higgs coupling a. Higgs self- 
coupling is taken as A = 10" 2 . Time is given in program units, 
see 12311 . Preheating ends earlier for Higgs field (t pr < 100) than 
for inflaton (t pr < 500). For a = 10 — 9 one has the border case 
when the average momenta of the fields are less then the lattice 
ultraviolet cutoff. 

The lattice results, however, show that such exponen- 
tial energy transfer into the Higgs particles for a broad 
range of parameters terminates before any significant 
part of the inflaton zero mode energy is depleted. The 
reason for that is the large Higgs boson self-coupling 
A ~ 0.1 which makes the re-scattering processes be- 
come important quite early. Unless the Higgs-inflaton 
coupling a is fairly large the re-scatterings terminate 
the resonance when only a negligible part of the energy 
in the inflaton zero mode is depleted! 4 1 On Fig. [2] one 
can see how the amount of the transferred energy de- 
pends on the value of the Higgs self coupling A which 
we allowed to vary to small values just to demonstrate 
the importance of the re-scattering processes. 

On Fig.[3]one can see the dependence of the total en- 
ergy transferred into the Higgs field as a function of the 
inflaton-Higgs coupling a. One can draw the conclusion 
that parametric resonance effects only become impor- 
tant at a ~ 10 -7 , which is too large a value. Thus, the 
reheating process proceeds by means of the simple de- 
cay of the inflaton (generated abundantly by parametric 
resonance) into the Higgs particle. This process will be 
analysed analytically in the next subsection, where we 
will advocate that this perturbative inflaton decay really 
reheats the Universe at lower values of the parameter a. 

2.1. Light inflaton case (mi < 2mn) 

While the parametric resonance regime for the Higgs 
is terminated quite early, the fluctuations of the inflaton 
field continue to grow exponentially. Since the amount 



4 For a potential without the inflaton mass term in a different part 
of the parameter space similar claims were made in fH . 
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Fig. 2. Energy transfer dependence on A (here 
a = /3 = 2.6 X 10~ 13 ). Values are taken at late time t pr = 10 3 . 
LatticeEasy parameters are as in Fig. [T] 
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Fig. 3. Energy transfer dependence on inflaton-Higgs coupling a. 
Values are taken at late time t E 
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show respectively the dependence of 

gives pff, ~ p x at a ~ 3 X 10~ 8 . LatticeEasy parameters are as in 
Fig. (TJ. For a physical value of the Higgs self-coupling (A ~ 0.1 
the energies become comparable even closer to a = 10 -7 . 



of the energy transferred into the Higgs field is practi- 
cally negligible the dynamics of the inflaton field is very 
close to that of the pure quartic inflaton model which 
was analyzed numerically in 12 ll 12211 . In brief, the in- 
flaton zero mode keeps driving the exponential grows 
of the nonzero modes until roughly half of its energy is 
transferred into the inflaton particles. After that the re- 
scattering processes become important, slowly moving 
the inflaton particle distribution to thermal equilibrium. 
At some moment the scattering process 21 — > 2H be- 
comes important and the Higgs particle (together with 
all other SM particles) is generated and the standard 
thermal history of the Universe takes over. The easi- 
est way to estimate the equilibration temperature of this 
process is to compare the mean free path na^i^H ~ 
-, where p avg is the average inflaton momentum, 



3 



with the Hubble expansion rateO H = ^r^ 1 ^- 
For the thermal distribution of the infiaton particles this 
leads to the estimate 



(8) 



However, the distribution of the infiaton excitations may 
be. generally, rather far from thermal equilibrium J21 




2211 . Evolution of the occupation numbers of the infiaton 
modes was found to be self similar in J21, 22] 



n(k, t) 



^o(fcr-f); 



(9) 



where r is the conformal time, k is the comoving mo- 
mentum, and p = 1/5 for three particle interactions and 
1/7 for four particle interactions, q ~ Ap. The only rele- 
vant for us property of the function no(fcr~ p ) is that the 
average momentum in at the beginning of reheating 
after inflation is l ^ i mY>\. Thus, the average momen- 
tum at later times is smaller than expected from the total 
energy density, p avg /T ~ (m P1 /T)P/?( 1+ P'/ 4 , where 
T ~ p 1 / 4 is now not a real temperature, but rather a pa- 
rameter defining the energy densitvPH (cf. equilibration 
time description in yj, 



2211 '). This enhances the 21 — > 
2H cross section together with the I number density, 
increasing the estimate (© by a factor (T/p avg ) 3 . This 
leads to the increase of the equilibration temperature by 
a factor 10 5 for four particle interaction, p = 1/7, and 
by a factor 10 2 for three particle interaction, p = 1/5. 
Exact calculation of the equilibration temperature re- 
quires extensive numerical study, but, in any case, the 
expression dS) should be considered as the lower bound, 
while 10 5 Tr is the upper (most conservative) bound. 

Requiring that Tr > 150 GeV we can obtain the 
lower bound on a 

a > 7.3 x 10" 8 , (10) 

for the thermal estimate dS) and 

a>7xl0~ 10 , (11) 

for the most conservative estimate of non-thermal dis- 
tribution of the inflatonFl 

While the bound (TTOb roughly coincides with the one 
at which the energy transfer to the Higgs field becomes 
effective enough to significantly deplete the zero mode 
of the infiaton (see Fig. |3]l while the value given by 
(fTTT i is about two orders of magnitude smaller. We can, 



5 topi = 2.44 X 10 18 GeV is the reduced Planck mass. 

6 After thermalization into the SM particles T transforms into the 
real temperature, up to the change of the number of d.o.f. 

7 Strictly speaking, one should also check if there is any kinematical 
suppression of the process. This may lead to 0(1) corrections and 
is, in fact, beyond the precision of present estimates. 



therefore, conclude that the upper bound on the infiaton 
mass is given by 



to/ < (0.14-5-1.40) 



/I 



150 GeV/ V 1.3 x 1Q- 13 



GeV , 
(12) 

where the range corresponds to the thermal or the most 
conservative non-thermal estimates. 



2.2. Heavy infiaton case (mj > 2m,u) 

In this case the infiaton mass allows for the direct 
decay of the infiaton into two Higgs particles. The cor- 
responding decay rate is given by 



2H) = \ 



2n 2 f3 



(3 mjj 
8ir m 3 j 



^-4 • (13) 



Comparing this rate with the Hubble parameter and 
requiring again for the reheating temperature Tr > 
150 GeV we get 



mi < 440 ^ 



4/8 



/3 



1/3 



GeV . 



-150GeV7 Vi-SxlO- 13 , 

(14) 

Of course, in the case a < (3/8 the generation of the 
cosmological perturbations is different from the case 
of pure quartic inflation. The Higgs field becomes rel- 
atively light and the parameter space of the model is 
modified. In particular, isocurvature fluctuations which 
one would generic ally expect in the two-field model 
have to be somehow suppressed. This will put the re- 
striction on the allowed values (a, f3). The analysis of 
this parameter space is very involved. One can expect, 
for example, that the parameter (3 can differ from its 
value in the case of pure quartic inflation. That is one 
of the reasons why the parametric dependence on j3 is 
kept in (flUiFI 

3. WMAP constraints and non-minimal coupling 

Finally let us discuss the constraints on the model 
from the WMAP data 111 611 . As was already mentioned 
in the inflationary regime the model is indistinguishable 
from the pure quartic potential inflation. One should 
then confront the fact that the amplitude of the tensor 
perturbations is too large. One possible resolution of 
this problem is to assume that the infiaton \ nas non- 
minimal coupling to gravity lfl7ll . We will repeat here 



8 Note, however, that the dependence of the bound 1141 on /3 is 
rather mild. 
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Fig. 4. The dependence of the quartic coupling j3 on the non-minimal 
coupling parameter £. 



the estimates following closely J25 . S0. We will 
take the following action as an example 



S = / d i xy / —g 



R 



■xd X ) 2 + m 2 -v(x,<$>) 



(15) 



where m ~ mpi. Even if the coupling £ is zero at 
a tree level one can expect that it will be generated 
via radiative corrections. As it will be discussed below 
even for small values of £ the coupling /3 will deviate 
from the one, obtained from the COBE normalization 
in the absence of the non-minimal coupling f3\^ = o ~ 
1.3 x 1CT 13 . 

The bound on the tensor-to-scalar ratio comes from 
the perturbations generated at N ~ 62 e-foldings (see, 
e.g. 1 150 before the end of inflation. In that regime the 
Higgs part of the model is not important and can be 
dropped to simplify the discussion. The inflaton part of 
(T5[ as it appears in Jordan frame by means of the con- 
formal transformation can be rewritten as (hat denotes 
transformed quantities) 



5.7 = 



d A x 



'PI 



R 



where 



g^v = Q g^ v , ft 2 ~ H T > 



'pi 



and the new field \ lS defined as 



pi 



dx V n 4 

The new potential is given by 



(16) 
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(18) 
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Fig. 5. Bounds on the inflation mass for successful reheating. Al- 
lowed regions for the case of II — » HH scattering (lower re- 
gion, light inflaton) and inflaton decay (upper region, heavy infla- 
ton). Higgs mass is chosen mu = 150 GcV and Tr > 150 GeV. 
Bounds are: I (inflaton decay), II (mj > 2mjj), III (2-2 scattering, 
non-thermal / distribution), IV (2-2 scattering, thermal / distribu- 
tion), V (a < 3 x 10~ 7 , smallness of radiative corrections). 



u(x) 



40(x)' 



-xHx) 



(19) 



We assume that £Xc/ to f>i ~ 1> where Xe is the value of 
the inflaton field at the end of inflation, so the contribu- 
tion to the effective Plank mass vanishes after the infla- 
tionary period. In that case the dynamics of the model 
with the action (T5[ after inflation is not different from 
that of the i/MSM model with the potential ©. This 
sug gest ion corresponds to £ < 0.1, see (fJTJ. Following 
1 171 |25J one can find that the first slow -roll parameter 
e is given by 



8mpi 



X 2 (m| 1 + £x 2 (l + 6£)) 



(20) 



Slow-roll ends when e = 1. From that one can find that 



1 



m pl 2(1 + 6£) 



^192f 2 + 32£ + 1 - 1 



(21) 



+ 0^), (£«!)■ 



The number of e-foldings from the moment when the 
inflaton field has the value xn till the end of inflation 
is given by 



TV = 



1 



'pi 



u 



(dU/dx) \d X 



f) dx 



(22) 



Xe 



(1 + 6£) - 6 In 



'pi 



'PI V '"PI + £Xe 

Since £ -C 1 one can find that with a good accu- 
racy xn ~ 2y 2 ^+(\P tnpi. The tensor-to-scalar ratio 
is given by ltl6ll 



5 



16e 



128m^ 



16(1 + 60 



1 



1) 



(23) 



(JV+l)(8£(iV- 

One can see lfl7ll that roughly in the interval £ = 
0.001 -e- 0.1 this ratio satisfies the WMAP constraints. 
The value of the inflaton self-coupling as a function 
of £ can be found from the COBE normalization 
U{xn)/^{xn) = (0.027mpi) 4 . The corresponding be- 
havior is shown in Fig. |4] This introduces slight growth 
of (3 with £, and thus increases all bounds simultane- 
ously, which is demonstrated in Fig. [5] 

4. Conclusions 

In Fig. [5] we combined the bounds on the inflaton 
mass we have found so far. We can conclude, therefore, 
that the mass of the inflaton in the z^MSM inflation |0] 
should be roughly in the range 

0.1 GeV < m/ < 10 GeV (24) 

in the case when it is light and in the range 

300 GeV < m, < 1000 GeV (25) 

in the case when the inflaton-Higgs coupling is very 
small. 

These bounds could be evaded in models with arbi- 
trary scalar field potentials, but the fact of the strong 
lower bound from reheating on the coupling between 
the inflaton and the Higgs should remain rather univer- 
sal. 

Values of £ larger then . 1 (and, therefore larger lower 
and upper bounds on the inflaton mass) are also allowed 
as well. However, since the dynamics of the model at 
preheating may be strongly modified from the one we 
have studied in this Letter it is hard for us to make 
any statements in that case, and we leave this for future 
analysis. 
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